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Abstract 

We suggest that recent neutrino puzzles that are the solar and atmospheric 
neutrino deficits as well as the possible neutrino oscillations reported by the LSND 
experiment and the possibility of massive neutrinos providing the hot component of 
the cosmological dark matter, can all be naturally explained by assuming existence 
of a mirror world described by an "electroweak" gauge symmetry [SU(2) x U(l)]' , 
with the breaking scale larger by about factor of 30 than the scale of the standard 
SU (2) x U(l) model. An interesting aspect of this model is that the sterile neutrinos 
arise from the hidden mirror sector of the theory and thus their lightness is more 
natural than in the usual neutrino mass scenarios. The needed pattern of the 
neutrino mass matrix in this model is obtained by assuming a conserved ZKM-type 
global lepton number L = L e + — L T , which is violated by Planck scale effects. 
One implication of our proposal is that bulk of the dark matter in the universe is 
a warm dark matter consisting of few KeV mass particles rather than the 100 GeV 
range particles of the currently popular cold dark matter scenarios. 



The present situation in the neutrino physics is rather intriguing but cautiously con- 
troversial. On the one hand, the direct measurements show no evidence for any of the 
neutrinos to be massive, providing only the upper bounds m Ue < 4.5 eV (< 0.7eV [2/3q„]) 
m v < 160 keV and m Ur < 20 MeV. The neutrinos could therefore be massless as far as 
these experiments are concerned. As is well known the standard model (SM) renders 
its left-handed (LH) neutrinos to be massless, since in absence of the right-handed (RH) 
states the lepton number conservation arises as an accidental global symmetry of the 
theory due to the joint requirement of gauge invariance and renormalizability. There is 
of course one conceivable scenario which could lead to neutrino masses within the min- 
imal SM: one could imagine that the minimal SM (or minimal SU(5) GUT) is a true 
theory up to the Planck scale, however the lepton number being a global symmetry, its 
conservation need not be respected by nonperturbative gravitational effects p], 0]. In 
that case, one can have nonrenormalizable operators of the type jj—l^jHH in the theory, 
where £1,2,3 are the lepton doublets and H is the Higgs doublet (In fact, this is equiva- 
lent to the seesaw mechanism || with the RH neutrino states having ~ Mpi Majorana 
masses.) These can then induce small neutrino majorana masses, at most of the order of 
m= (H) 2 /Mpi = 3 • 10~ 6 eV 0, which value can be naturally considered as a neutrino 
mass unit in the SM. 

On the other hand, there are indirect "positive" signals for neutrino masses and mixing 
accumulating during past years. In particular, if any of the following hints will prove to be 
true, this would point towards neutrino masses much larger than m. This hints include: 

(a) The solar neutrino problem (SNP). The solar neutrino experiments [|J] indicate 
a deficiency of solar v e fluxes which cannot be explained by astrophysical reasons [gj. 
This then implies that the discrepancy between theoretical expectations and experimental 
observations is due to new neutrino properties, the most plausible one being the oscillation 
of v e into another neutrino v x . Interestingly, the long wavelength "just-so" oscillation with 

rt O 2 

5m ex r^m and large mixing angle, just as the Planck scale effects predict, can provide 
a SNP solution ||. However, the MSW oscillation appears to be the most popular 
and natural solution to SNP, based on v e — v x conversion in solar medium. The required 
parameter range corresponds to 5m 2 ex ~ 10~ 5 eV 2 and, in the preferable "small mixing 
angle" scenario, sin 2 28 ex ~ 10~ 3 — 10~ 2 . 

(b) The atmospheric neutrino problem (ANP). There is an evidence for a significant 
depletion of the atmospheric z/ M flux by almost a factor of 2 ||. This points again to 
neutrino oscillations, with a parameter range sin 2 29^ x ~ 1 and 5m 2 ~ 10~ 2 eV 2 . 

(c) Dark matter problem. The COBE measurements of the cosmic microwave back- 
ground anisotropy suggests that cosmological dark matter consists of two components: 
cold dark matter (CDM) being a dominant component (^cdm — 0.7) and hot dark mat- 
ter (HDM) being a smaller admixture (^hdm — 0.25) || [TOfl . The latter role can be 
naturally played by neutrinos with mass of few eV's. As for the CDM candidates, several 
possibilities can be envisaged. For example, in supersymmetric theories with conserved 
R parity it can be provided by lightest neutralino. However, it can be of interest to 
think of CDM as also consisting of neutrinos, this time heavier (keV range), and with 
correspondingly small concentration, so called warm dark matter (WDM) |jTT]| . 

(d) LSND result: Direct evidence of — v e oscillation from the recent Los Alamos 
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experiment [0, with 5m 2 efl > 0.3 eV 2 and sin 2 26> eA1 = 1(T 3 - 10~ 2 . 

Thus, these hints if proved correct would point to the neutrino masses much larger than 
m, which in turn would imply nontrivial physics beyond the Standard Model at a scale 
below Mp\. In particular, this would imply that the neutrino masses are induced from 
the effective operators of the type ^kljHH with the cutoff scale A <C Mpi. Since this 
operators with A -C Mpi cannot be seriously considered as truly nonrenormalizable, they 
have to be effectively generated via some renormalizable couplings of certain additional 
particles - fermions or scalars, and A is related to their mass scale. For example, in 
the context of seesaw mechanism || these are the RH neutrinos. Other neutrino mass 
mechanisms as are the Zee model, the heavy triplet mechanism etc. (see e.g. |HJ and 
references therein), also effectively reduce to the operators of these type. 

A pessioptimisticQ approach to the situation would be to assume that all these indica- 
tions for neutrino masses will indeed be borne out by future experiments and one should 
therefore study its implications for the texture of neutrino masses. Such a study has been 
carried out in the past two years [14] and has led to the amazing result that only one pos- 
sible neutrino mass texture is compatible with all the above mentioned data. It requires 
an extra light, sterile neutrino v s beyond the three known neutrinos (i^ e , and v T \ In 
this scenario the SNP is explained as a consequence of v e -v s MSW oscillation while the 
atmospheric neutrino puzzle is explained by — v T oscillation with m u ~ m VT ~ 2.4 eV. 
One assumes in this scenario that m v <C rn v . The and v T provide the cosmological 
HDM and it also explains the LSND results. The detailed mass matrix for the neutrinos 
in this case is essentially in a block matrix form with each block being 2x2 with a small 
mixing between the v s and to account for the LSND result. 

Indeed, the small mixing angle MSW oscillation v e — v s is the only place where the 
extra sterile neutrino could show up. Obviously, sterile neutrino cannot constitute HDM, 
neither it can be applied to ANP, due to primordial nucleosynthesis bound on the effective 
number of neutrino species at t ~ 1 s: N u < 3.1 (|IJ|. In particular, this bound implies that 
the oscillation of active neutrinos into the sterile one should obey the limit Sm 2 sin 2 29 < 



1.6 ■ 10 6 eV 2 [19 , which excludes the — v s oscillation as a candidate for the ANP 



solution. The same bound excludes the possibility for the SNP solution through the 
MSW oscillation v e — v s in the large mixing angle regime. 

Concerning the sterile neutrinos, an immediate question that arises is 'where do they 
all come from, and where do they all belong': namely, why they can be so light if their 
masses are allowed by the gauge symmetry. In the SM a role of sterile neutrinos can be 
played by the RH components of the usual ones. However, in this case it is difficult to 
realize why they in combination with LH neutrinos do not form the Dirac particles as 
heavy as the charged fermions, or why they do not have large Majorana masses in the 
spirit of the seesaw scenario. In order to render sterile neutrinos massless, some ad hoc 



global symmetries should be introduced WE, [L6| . The recent proposal |T7| that the sterile 



neutrino can be related to axino or majorino like fields in the context of supersymmetric 
theory is also based on ad hoc constraints. 

Here we suggest that the sterile neutrinos are in fact the neutrinos of a mirror world 
which is the mirror duplicate of our visible world except that its " electroweak" scale 



^Here our points: pessimistic (Z.G.B) and optimistic (R.N.M) did not agree. 
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v' is larger by factor of ( than the standard electroweak scale v. It interacts with the 
visible world only through gravity and possibly through some superheavy singlet scalars. 
Thus, neutrinos of the mirror standard model v' e n T should be light by the same reason 
as ordinary ones v' e T - their direct mass terms are forbidden by the mirror SM gauge 
symmetry. Then the SNP can be explained through the oscillation of v e into its mirror 
partner v' e by means of Planck scale effects, which naturally provide both the values of 
5m 2 and sin 2 29 in the range needed for the small mixing angle M SW solution. The 
question arises, what is a possible role played by the two remaining mirror neutrinos. Our 
answer is that in the framework presented below they have masses in the KeV range and 
constitute the warm dark matter of the universe. We also address their implications for 
nucleosynthesis. 

While the SNP solution can be completely understood by means of Planck scale in- 
duced MSW oscillation v e — u' e , understanding of other puzzles asks for some non- minimal 
mechanisms generating masses of the v^ T states. We suppose that the neutrino sector 
of the model obeys the separate lepton number conservation in both sectors, and assume 
that these lepton numbers are broken in two stages: (A) the dominant entries (~ eV) in 
the neutrino mass matrix has origin in some intermediate scale physics which respects a 
global ZKM-type lepton number (L = L e + L M — L T in our model); this fixes the skeleton 
of neutrino mass matrix with specific texture needed for reconciling the LSND oscillation 
and the HDM; (B) the conservation of this global number is violated by Planck scale 
effects, which can then explain SNP and ANP. In particular, we show that if the scale of 
the mirror world is some 30 times the electroweak scale, then the Planck scale induced 
terms reproduce the orders of magnitudes of the parameters required by the small mixing 
angle MSW scenario. 

In the spirit of our hypothesis, there will be a similar pattern for the Majorana mass 
matrix in the mirror sector except that the entries should be scaled up by a factor of £ 2 . 
As far as ordinary neutrinos v^^ with ~eV mass play a HDM role, their mirror partners 
v 1 T being ( 2 times heavier, for £ ~ 30 fall in keV range and thus can constitute WDM of 
the universe. 



Without the mirror world 



Apparently, three known neutrinos are not enough to reconcile all these hints |14| . 
In fact, the key difficulty is related to SNP. The points (b), (c) and (d) can be easily 
reconciled e.g. by assuming that in the basis of the flavour eigenstates v e ,fi,T the neutrino 
Majorana mass matrix has the texture obeying the L = L e + L M — L T conservation: 




(1) 



This matrix has one massless (Vi) and two massive (^2,3) degenerated eigenstates, with 
mass m2,3 = m = (a 2 + 6 2 ) 1 / 2 . The latter can play a role of HDM provided that m is 
of few eV's (e.g. m ~ 2.4 eV, according to ref. 
mixed by the angle 9 e ^ (tan^ = a/b), and 5m 
LSND oscillation if sin 2 20 e/i ~ 2 • 10~ 3 @. 
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On the other hand, v e and are 
m 2 ~ 6eV 2 which can explain the 
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Small explicit violation of L would induce nonzero entries in the matrix ([I]) and thus 
triger the — v T oscillation. For example, if the small nonzero entry e <C m appears in 
its (3,3) element, then we have 5m 2 T « 2em and almost maximal mixing, sin 2 29^ T « 1. 
Thus, APN can be solved provided that 8rn 2 ~ 10~ 2 eV 2 , which for m ~ 2.4 eV in turn 
implies e ~ 2 ■ lCT 3 eV. 

Thus, only the SNP remains unresolved. 



Connecting the Two Worlds 

Imagine that besides the world of the usual particles of the standard electroweak model 
G = SU(2) x i7(l), there exist an analogous world of particles belonging to the mirror 
gauge group G' = [577(2) x {7(1)]'. In particular, the ordinary and mirror lepton states 
and Higgs doublets transform as 

G: k= ( ^ ~ (1,-1; 0,0), Ef ~ (0,2; 0,0), if~(i,l;0,0) 
G': l\= ^ U j_ J ~ (0,0; 1,-1), Ef ~ (0, 0; 0, 2), ff' ~ (0, 0; |, 1) (2) 

where the weak isospins I, I' and hypercharges Y, Y' are shown explicitly, and index 
% = 1, 2, 3 denotes the electron, muon and taon families respectively. One can impose also 
the invariance under discrete transformation P(G <-> G') simoultaneously interchanging 
all corresponding particles of the ordinary and mirror worlds.0 We assume that P is 
spontaneously broken, for example by the VEV of some P-odd singlet scalar 77 coupled to 
H and H' As a result, H and H' can have different nonzero VEVs: v' ^> v — 174 GeV. 
Hence, the mirror world is completely analogous to ours, but with all particle masses just 
scaled by the factor ( = v'/v. In particular, for the gauge bosons we have Mw',z> = (M wz 
and, as far as the charged leptons get masses through the P-invariant Yukawa couplings 
giliEfH + gil^EfH', we also obtain m' = (m e ^ )T . The same applies to the ordinary 
and mirror quark masses. Photons remain massless in both of these worlds |ZD| . 



Let us suppose that the two worlds comunicate only through gravity and possibly also 
via some superheavy gauge singlet matter, like the P-odd scalar rj. We do not consider the 
possibility of the mixed representations. Concerning thermodynamics of the two worlds 
in the Early Universe, we assume that at the inflationary reheating temperatures they 
are already decoupled from each other. If the inflaton couplings violate P-invariance, 
then one can imagine the situation when the visible and mirror particles are "reheated" 
with different rates, so that after inflation the effective temperatures of the ordinary and 
mirror thermal bathes are different. As a result, present cosmo logical abundance of the 
mirror particles (including lightest "charged" states e', u', d! as well as mirror neutrinos 
v' e T and photons 7') can be less than that of their visible partners. The contribution of 

2 ) All fermion states are given in the LH basis. Quarks can be included in a strightforward way: the 
ordinary quarks transform as triplets of SU(3) C while the mirror ones are triplets of the mirror gauge 
group SU(3)' C . Notice that such a theory can be consistently extended to the SU(5) x SU(5)' GUT, also 
in the supersymmetric version. 

3 ' In fact, the ordinary (LH) fermions can be interchanged with the mirror fermions of the same helicity, 
as well as with their conjugated (RH) states. In the latter case, P formally plays a role of parity. 
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v' and 7' to the universe expansion rate at the nucleosynthesis epoch t ~ 1 should be very 



small: N u < 3.1 which can be translated into the upper limit on today's abundance 
r = n v >/n v < 0.04 - 0.02 for ( varied from 10 to 100.0 

With the minimal particle content (0), the neutrinos of both sectors stay massless 
unless one appeals to the gravity induced Planck scale effects which explicitly violate the 
global lepton number, and also can mix the neutrino states of the two worlds ||. The 
relevant higher order operators are (C matrix is omitted): 



^ikm.H) + j^mw.H') + ^ium-H') + h. c . ( 3 ) 

with the constants a, (3 ~ 1. More in general, an order of magnitude less values of these 
constants should not come as a surprise, while it is also possible (e.g. in the context of 
string philosophy) that the actual cutoff scale of the operators (|) is order of magnitude 
less than Mp\. Hence one can let these constants to range say from 0.1 to 10. Obviously, 
these operators generate the neutrino mass terms in both worlds, as well as their mixing 
term. Indeed, after substituting the VEVs of the H, H' we obtain 

m *» = a TF~ = am i m K = a TT~ = C am, m Ve< = fJ—— = (f3 m (4) 
Mpi Mpi Mpi 

where m= 3 • 10~ 6 eV. As a result, v e — v' e oscillation emerges which can explain SNP: 

Sm 2 = a 2 (C/30) 4 x 8 • 10~ 6 eV 2 , sin 2 26 = ^(30/C) 2 x 4.5 ■ 10" 3 (5) 

a 2 

For a, (3 ~ 1, the choice ( ~ 30 fits the parameter range in the regime of "small mixing 
angle" MSW solution 0. More in general, taking into account the solar model uncertain- 
ties in predicting the boron and beryllium neutrino fluxes, the relevant parameter range 
can vary within sin 2 26 = 6 • 10~ 4 - 2 • 10^ 2 and 5m 2 = (4 - 10) • 10" 6 eV 2 [0]. Then, by 
assuming a ~ /3, the above range for sin 2 26 can be recovered for ( = 10 — 100, while the 
mass difference is correlated with sin 2 26 as 

5m 2 „ a 2 ( t^Ll) 2 x 6 . 10-6 eV 2 ( 6 ) 

V sin 2 26 J K ' 

which is compatible with the MSW range for a proper a in the interval 0.1 — 10. 

The explanation of the other neutrino puzzles needs larger masses, which cannot be 
born by the Planck scale effects. Thus, some non-minimal neutrino mass mechanism 
should be incorporated. As we stated earlier these masses in each sector can arise only 
through the effective operators respectively bilinear in H and H' . Such operators can 



4 ) These estimates depend on the abundance of v' relative to at t ~ 1 s. The neutrino decoupling 
temperatures in the mirror and ordinary worlds are related as T' d /Tei ~ C 4 ^ 3 : while the fermion masses 
are scaled by factor £. E.g. for ( ~ 10, v' decouples from 7' after the QCD' phase transition (due to P 
invariance, the 'mirror' confinement scale cannot be substantially larger than that of the visible world: 
Aq CD < lGeV). Then for their effective temperatures at t ~ 1 s we obtain T v i /Ty — (4/11) 1 / 3 . On the 
contrary, for ( ~ 100 the decoupling temperature T' D > 1 GeV and thus the contribution of light mirror 
quarks u',d' has also to be taken into account, which implies T v * /Ty = (4/53) 1 / 3 . 
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be effectively induced by several mechanisms and we will discuss one such mechanism in 
brief in the next section. But on general grounds, one can write 

^{kH^H) + ^{l' i H'){l' j H') (7) 

where are some 0(1) 'Yukawa' constants with the L-conserving texture (|1]), and A ~ A' 
are certain regulator scales.Q Then the V\ ~ v e state and its mirror partner u[ rs v' e are 
rendered massless so that the SNP solution can be ascribed to the Planck scale effects in 
a manner demonstrated above, provided that ( ~ 30 or so. As for the degenerated states 
z/ 2; 3 and ^2 3 , their masses are scaled as vnl ~ ( 2 m. Therefore, for m in the eV range as 
it is needed for explaining the origin of HDM, u' 2 3 emerge in the keV range and thus can 
constitute the WDM. 

Let us remark that the ANP remains unresolved unless one introduces the explicit 
mechanism for L breaking providing nonzero entries e ~ 10~ 3 eV in the mass matrix ([J). 
The Planck scale operators (H) can provide at most ~m contribution which is far below 
the needed value. However, in the next section we demonstrate a model in which APN 
solution can be also related to the Planck scale physics. 



Origin of the Dominant Neutrino Masses 

Let us briefly discuss how the dominant neutrino mass matrix for each sector given 
in Eq. ([[]) can emerge. As we already stated, the neutrino masses in each sector can 
arise only from the operators respectively bilinear in H and H' . These operators can 
be effectively induced by several mechanisms. The model we consider is a variant of 
the so called /z-model P2fl . Let us assume that theory obeys conservation of individual 



global lepton numbers Lj = L e ^ T prescribed to lepton states in (^j in obvious way, which 
is spontaneously violated by the VEVs of the gauge singlet scalars with mixed lepton 



numbers: $ eAi , <£> er and ^i-, in the sprit of refs. |23|, [ilj. In particular, we assume that 
only ($ eT , $ M t) ~ V are nonzero, while ($ eAt ) = 0. In this way, the ZKM-type lepton 
number L = L e + — L T remains conserved. We also introduce the additional neutral 
fermions A^(— 1), Si(+1) and iV|(— 1), S'^+l), where brackets show the individual lepton 
charges L«. Then the G x G' x P invariant coupligs, also respecting the conservation of 
lepton numbers, are the following:0 

Cvuk = MiNiH + I'MH') + M,(A^ + N'iSfi + A<% (S^ + S'^) (8) 

For simplicity we assume that the Dirac mass terms are closely degenerated: Mi.2,3 ~ 
M 3> v', while the Yukawa couplings fi have approximately the same pattern as that of 
the charged fermions, say /3 ~ 1, f 2 ~ 3 • 1CT 2 and fi ~ 10 -3 . 

5 - ) We will assume these scales to be independent of the VEVs of H and H' . For simplicity, let us take 
A' ~ A, which seems especially natural if A ^ v',v. Indeed, in order to induce ~eV entries in the mass 
matrix m„ (|l|), A should be of about 10 13 GeV. 

6 ' Notice however that not all terms allowed by the symmetry are introduced. Their absence can be 
ensured by imposing a conservation of the individual global lepton numbers L e , L^, L T and L' el L' L' T 
separately in each sector. Alternatively one can think of the N, S and N', S' states as being the SU(2) 
and SU(2)' isotriplets with zero hypercharges. For example, in the 51/(5) x SU(5)' extension of our 
model these can be originated from the heavy fermionic 24-plets in each sector. 
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Then the total mass matrices of the neutral states in each sector have the form 
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(9) 



where the Majorana mass matrix /i has a texture conserving L 

/ A 



B 
\AB0 



A~ B ~\V 



(10) 



Let us now take into account also the effects of the possible higher order operators 
explicitly violating the lepton number conservation. Among the plethora of possible 
Planck scale operators relevant are only the ones given by eq. (|j) and 



M, 



(SiSj 



in 



where the <fi 2 stands for gauge invariant bilinears of the scalar fields involved and obviously 
it is dominated by the scalars $jj having the largest VEVs in our model. Thus, the mass 
matrix of the S states becomes /t + 7 e, where e= V 2 /Mpi. As far as the latter operator 
violates L too, the matrix 7 in general has no zeroes. 

Then, after decoupling the heavy states the mass matrix of the light neutrinos becomes 



m,. 



s(/t + 7 e)s+ m a ( m {3 

Cm/3 T C 2 S(A + 7£)s + C 2 



m a 



§ = diag(si,s 2 ,s 3 ) (12) 



where Sj = fiv/Mi. Clearly, in the absence of the Planck scale induced corrections the 
mass matrix m v of the active neutrinos has the L conserving texture (P with a = S1S3A 
and b = S2S3B, while the mass matrix of their mirror partners is just scaled by factor of 
( 2 rh u . Planck scale induced corrections violate L, giving rise to small nonzero 



C 

entries in neutrino mass matrix. Let us consider in more details how this model solves 
the present neutrino puzzles: 

HDM+WDM. In order to represent the cosmological HDM component, the mass m 
of the degenerated states u 2 ,3 should be in the eV range, say m = 2.4 eV [Kj. Then for 
C ~ 30 masses of their mirror partners m' = ( 2 m are in the keV range and thus can 
constitute cosmological WDM provided that their present abundance relative to active 
neutrinos has a proper value. This can be expressed as 



2m = n HDM • 92/i 2 eV, 2rm = tt WDM ■ 92h 2 eV 



(13) 



where h is the Hubble constant in units 100 Km s _1 Mpc -1 . Cosmological density of 
the visible baryonic matter corresponds to Qb — 0.05. We also recall that in our model 
due to different inflationary reheating, VLb> <Qb unless baryon asymmetry in the mirror 
world is much larger than that of normal baryons. Thus by taking rather conservatively 
Qb + Q' B < 0.1, the remaining cosmological density can be shared between HDM and 
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WDM components as say Qhdm — 0.2 and Ohdm — 0.7. Then for the needed abundance 
of the mirror neutrinos relative to the active ones we obtain r =~ 3.5 • (~ 2 . For ( > 10 
this is obviously consistent with the upper bound from the nucleosynthesis constraint 
r < 0.04- 0.02. 



LSND: Oscillation v e — occurs with = m 2 ~ 6 eV 2 and sin 6*^ = f\Aj f^B. 

Thus for our pattern of the Yukawa constants /; we obtain sin 26 efl ~ 10~ 3 , in the range 
needed for for the LSND oscillation, 

ANP: The Planckian terms (|TT|) induce nonzero entries in the matrix (x of the order of 
£— V 2 /Mpi. This in turn induces in rh v nonzero (3,3) entry e = S37 e which removes the 
b> 2 — v<i degeneracy and thus trigers — v T oscillation with sin 29 « 1 and Sm 2 ^ = 2em. 
Since m = S2S3XV ~ 2.4 eV, then the atmospheric neutrino oscillation range requires 
e ~ 2 ■ 10~ 3 eV. Thus, confronting the values of m and e, we obtain the following estimates 
for the lepton number breaking scale V and the Dirac mass M: 

V = ^M Pl ~ * • 3 • 10- GeV, M = • ~^=M Pl ~ ' 1Ql3 GeV ( 14 ) 

/37m 7 hT /2 m s/eMpx 7 1/2 

so that both scales are large for reasonable values of A and 7. 

SNP: Obviously, the operators (|Tl~|) negligibly contribute to the electron neutrino mass. 
Indeed, their contribution to the (1,1) entry in the matrix m v is (fi/ fz) 2 £ <^?ri. Thus, 
for the mass terms of v e and v' e the relevant contributions come dominantly from the 
operators @, leading to MSW oscillation with the parameter range given in (|5]). 

The following remark is also in order. Once the individual global lepton numbers 
L e ,p,T are broken spontaneously down to L = L e + L M — L T , there should exist flavons 



23| , [16fl , multiflavour analogues of the singlet majoron [gj]. As is known, the Planck scale 
effects applied to majoron-flavons, induce their masses ~ (V / M) l l 2 V which overdose the 
universe if the massive majorons present in the today's universe f25|. The estimates ( |i~4|) 
demonstrate, however, that the scales V and M are enough large, and appearance of the 
majoron after the inflationary reheating can be avoided. 

In conclusion, we have shown that if the existence of the light sterile neutrino is to 
be taken seriously for solving the neutrino puzzles, then the its lightness is most easily 
understood if it is identified with the mirror neutrino. We have then outlined a two step 
mechanism whereby the desired texture needed to solve the neutrino puzzles in a unified 
manner can emerge. We have given also an explicit model that can lead to the basic 
texture. 



Acknowledgements. Z.G.B thanks V. Berezinsky, A. Dolgov, G. Fiorentini and Ursula 
Miscili for valuable discussions. The work of R.N.M is supported by a grant from the 
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Note added: When our work was finished, we encountered the paper by R. Foot and 
R. Volkas, preprint UM-P-95/49 (May 1995), in which the explanation of recent neutrino 
puzzles is attempted with the mirror model with exact parity. Our approach and results, 
however, are drastically different from those obtained in this paper. 
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